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^ ABSTRACT 

_^ 

^^ We report the detection of CoRoT-18b, a massive hot jupiter transiting in front of its host star with a period of 1.9000693 ± 

l/~J 0.0000028 days. This planet was discovered thanks to photometric data secured with the CoRoT satellite combined with 

T-H spectroscopic and photometric ground-based follow-up observations. The planet has a mass Mp= 3.47 ± 0.38 Mjup, a radius 

Rp = 1.31 ± O.lSRjup, and a density pp — 2.2 ± 0.8g/cm^. It orbits a G9V star with a mass M* — 0.95 ± O.ISM©, a radius 

I I Rt = 1.00 ± O.I3R0, and a rotation period Prot = 5.4 ± 0.4days. The age of the system remains uncertain, with stellar evolution 

Ph models pointing either to a few tens Ma or several Ga, while gyrochronology and lithium abundance point towards ages of a few 

PJ hundred Ma. This mismatch potentially points to a problem in our understanding of the evolution of young stars, with possibly 

significant implications for stellar physics and the interpretation of inferred sizes of exoplanets around young stars. We detected 

the Rossiter-McLaughlin anomaly in the CoRoT-18 system thanks to the spectroscopic observation of a transit. We measured the 

obliquity ^p = 20° ± 20° (sky-projected value A = —10° ± 20°), indicating that the planet orbits in the same way as the star is 

Q rotating and that this prograde orbit is nearly aligned with the stellar equator. 

\^ 
-i-^ Key words, stars: planetary systems - techniques: photometry - techniques: radial velocities - techniques: spectroscopic 

1. Introduction lows several thousand stars (V^ = 12 — 16) to be continu- 

^^ ously observed for up to 150 days with a high photometric 

> Out of the ^ 550 exoplanets known to date, more than 100 accuracy CoRoT is thus well adapted to detecting tran- 
CN transit their parent stars as seen from the Earth. This par- ^^^-^^ ^^^^^^^ ^-^^ g^^^U ^^^^-^ ^^^-^ ^^ CoRoT-7b (Leger 
^ ticular configuration allows numerous key studies, including ^^ ^i_ ^^ Q^gj^^ ^^ ^1 ^qq^ ^ ^^ ^^ ^^^^ ^^^i^^l periods, 
>0 accurate r adius m ass, and thus density measurements (see, ^^^^ ^ CoRoT-9b (Deeg et al. mM ■ It can also detect hot 

. e.g., Winn|2010a|for a review), atmospheric studies m ab- j^piters, such as CoRoT-18b. We report its discovery here. 

r- sorption through transits and memission through occulta- ^^ describe in Sect. El the CoRoT observations and 

O tions (e.g. Vidal-Madjar et al. [20031, Wheatley et al. [2^01), ^^e transit detection of the planetary candidate. Then, we 

Z2 dynamic analysesthrough possible timing variations (e.g. p^^^^^^ ^^ Sect. [31 the ground-based follow-up observations 

. . Holman et al.|MnD, or spm-orbit alignment measurements ^^^^ ^^^^ needed to establish the planetary nature of the 

> *^^S^° * "" Rossiter-McLaughlm effect (e.g. Bouchy et ^^^^^ detected by CoRoT and also to characterize this plan- 
al. mm- The power of these analyses incited numerous ^^^^y ^^^^^^^ The analysis of the whole dataset and the re- 
search surveys for transiting planets. Most of them were ^^j^g ^^^ presented in Sect. [1 before conclusion in Sect. § 

J3 discovered in the last five years, and the detection rate is i_i 1— 1 

still increasing. 

The CoRoT space mission {COnvection Rotation and 2. CoRoT observations and transit detection 
planetary Transits, Baglin et al. 2003^ was launched on 

2006 December 27. Based on a 27-cm telescope and a CoRoT-18 was one of 4161 target stars observed by CoRoT 

2.8° X 2.8°-field camera, it is designed to study asteroseis- ^om 2010 March 5 to 29 as part of SRa03, the third short 

mology and detect transiting exoplanets. The satellite al- run of the satellite in the Galactic anti-center direction. 

The coordinates, magnitudes and identifiers of CoRoT-18 

Send offprint requests to: G. Hebrard (hebrard@iap.fr) ^ various catalogs are given in Table [l] The finding chart 

* The CoRoT space mission, launched on 2006 December 27, i^ plotted in Fig. ^ Following the method described in 

has been developed and is operated by CNES, with the con- Gandolfi et al. (2008), the distance, d, and interstellar ex- 

tribution of Austria, Belgium, Brazil, ESA (RSSD and Science tinction. Ay, to CoRoT-18 were derived using the DENIS 

Programme), Germany and Spain. and 2MASS magnitudes reported in Table [ll and synthetic 
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Table 1. CoRoT-18 IDs, coordinates and magnitudes. 



CoRoT window ID 


SRa03_E2_1347 


CoRoT ID 


315211361 


USNO-Bl ID 


0899-0092144 


2MASS ID 


06324137-0001537 


CMC14 


063241.3-000153 


Coordinates 






RA (J2000) 


06'' 32™ 


41!36 


Dec (J2000) 


--00° 01 


53'.'71 


Magnitudes 






Filter 


Mag 


Error 


B2 (USNO-Bl) 


15.79 


- 


R2 (USNO-Bl) 


14.99 


- 


V 


15.00 


0.10 


R (CMC14) 


14.472 


0.048 


I (DENIS) 


14.051 


0.030 


J (2MASS) 


13.441 


0.024 


H (2MASS) 


13.080 


0.031 


K (2MASS) 


13.014 


0.030 



9^ 



^j 



Fig. 1. The sky area around CoRoT-18 on the POSS red 
image. The target is in the middle of the image, with the 
overplot of the CoRoT photometric aperture mask. 



colors from a model at mos phere with the same parameters 
as the star (see Sect. 4.2 and Table [3 below). We found 
Ay = 0.15 ± 0.15 mag and d = 870 ± 90 pc, as well as 
V = 15.00 ± 0.10 in the Johnson standard system. 

The cadence for this target was 32 seconds throughout 
the CoRoT observations, resulting in 65120 exposures span- 
ning the wavelength range 300-1000 nm. The data were pro- 
cessed by the CoRoT pipeline (Auvergne et al. 2009). The 
normalized white-light curve, obtained by summing the flux 
from the three channels and normalizing by the mean flux, 
is shown in the top panel of Fig. [2J It clearly shows 13 
planetary-like transit features, with a period of ^ 1.9 days 
and a depth of ~ 2 %. The phase- folded light curve is plot- 
ted in the upper panel of Fig. [S] 

To reduce computing time, the data was rebinned to 
512 seconds, which is the normal observing cadence for 
CoRoT exoplanet targets. This binned light curve is used 
throughout the rest of this paper. We checked that this 
binning does not significantly affect the parameter retrieval 
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Fig. 2. Top: CoRoT light curve in black with the 13 tran- 
sits shown in red. All epochs in this paper are given in 
Heliocentric Julian Date (HJDutc; Eastman et al. I2010p . 
The original light curve (shown in gray) contained a discon- 
tinuity at HJD -2450000 = 5281.62, which was corrected 
by subtracting 15 mmag from data taken after this date. 
A linear fit to the light curve has also been used to re- 
move any trend on timescales longer than the duration 
of the run. The light curve overplotted in blue is binned 
to one point per orbital period of the CoRoT satellite to 
make sure that no systematics from the rotation period can 
come through. Bottom: Lomb-Scargle periodogram of the 
out-of-transit light curve (black curve in top panel) as a 
function of relative improvement in x^ compared to a con- 
stant flux model. It shows the stellar rotation signature at 
5.4 ± 0.4 days (solid red line), as well as its first harmonic 
(dashed red line). Also shown are the Gaussian fit used to 
estimate the uncertainty in the stellar rotation period and 
the resulting l-a interval (dotted red lines). 



by using the formalism presented by Kipping (20 10b|) . The 
mean fiux was 57366.4 e~ per 32-second exposure, and the 
relative standard deviation of the binned light curve is 
6.9 X 10~^, a factor 7.2 above the photon noise. This factor 
reduces to 1.6 when the transits are removed and variations 
on timescales longer than a day are filte red o ut (black line 
in upper panel of Fig. [2J see below Sect. 4.1 1. 

CoRoT-18 was one of ten objects of interest identified 
soon after the end of SRa03 observations by the "alarm 
mode" pipeline (Surace et al.[2008), which removes the out- 
liers flagged by the main pipeline, detrends the light curves 
to remove long-term trends (instrumental and stellar) using 
a median filter, and searches for transits using an implemen- 
tation of the box least squares (BLS) algorithm of Kovacs, 
Zucker & Mazeh (2002). A number of tests were then car- 
ried out to check that the CoRoT data were compatible 
with a planetary origin for the transits: full transit fits to 
the white, red, green, and blue light curves using the transit 
formalism of Mandel & Agol ()2002|1 . check for differences 
in depth between odd- and even-numbered transits, search 
for possible companion occultation at the transit antiphase, 
and search for ellipsoidal variations. Because this candidate 
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Fig. 3. Phase-folded CoRoT (top panel) and Euler (bot- 
tom panel) transit light curves with the best fit and the 
residuals. CoRoT data cover 13 transits with a 512-second 
bin and the Euler data one transit with a ^ 3-minute bin. 



Both datasets are joint-fitted (see Sect. 4.1) 



passed all these tests, it was put forward for follow-up ob- 
servations with high priority. 

3. Ground-based follow-up observations 

3.1. "On-ofr photometry 

The point spread function of CoRoT that contains 50% of 
the fiux has an elongated shape of about 35" x 23". The pho- 
tometry is done through an aperture of that size. Owing to 
this poor spatial resolution, a deep stellar transit diluted in 
the flux of other source(s) also included in the large CoRoT 
aperture could mimic a shallow planetary transit. "On-off" 
photometry of the transit performed from the ground with 
a telescope allowing higher spatial resolution could identify 
contaminated eclipsing binaries (Deeg et al. I2009p . 

"On-off" photometric observations of CoRoT-18 were 
performed in November 2010 with the ESA Optical Ground 
Station (OGS) 1-m telescope, located at Izaha in Tenerife 
(Spain). Alternated short exposures of 30 seconds and long 
exposures of 100 seconds were taken, for a total dura- 



tion of 28-min observation on-transit and 28-min obser- 
vation off-transit. Seven extra sources are detected in a 
radius of 30" around the main, brighter target. We per- 
formed aperture photometry of the target and neighbor- 
ing stars. The main target shows a transit with a depth 
0.03±0.01 mag, in agreement with the transit detected with 
CoRoT. The seven other targets show stable fluxes, within 
precisions ranging from 0.01 to 0.39 mag, depending on the 
target and the exposure time. This "on-off" observations 
thus allowed us to exclude the detected transit signature 
caused by an eclipsing binary diluted in the CoRoT point 
spread function. 

3.2. Radial velocities 

We started the spectroscopic follow-up of CoRoT-18 in 
October 2010 with the SOPHIE spectrograph at the 1.93- 
m telescope of Haute-Provence Observatory, France. Three 
measurements performed in three successive nights near ex- 
treme phases (assuming a circular orbit) showed large radial 
velocity variations, in phase with the CoRoT ephemeris. 
The variation was on the order of 1 kms^^, indicating a 
companion with a mass around 3Mjup. Thus wc decided 
to pursue the spectroscopic observations with SOPHIE to 
strengthen the detection and to characterize the planetary 
system. We also used the HARPS spectrograph at the 3.6- 
m ESO telescope in La Silla, Chile, and the 2.56-m FIES 
spectrograph attached at the Nordic Optical Telescope in 
La Palma, Spain. Having three ground-based instruments 
at different longitudes was useful for reaching a good phase 
coverage for this system, which has an orbital period close 
to an integer number of terrestrial days. The observa- 
tions were conducted up to January 2011, in good enough 
weather conditions to allow satisfactory data to be secured 
in reasonable exposure times for this faint target. 

B othSO PHIE (Bouchy et al. [20091 and HARPS (Mayor 
et al. I2003P are cross-dispersed, environmentally stabilized 
echelle spectrographs dedicated to high-precision radial ve- 
locity measurements. SOPHIE data were acquired in High- 
Efficiency mode (resolution power R — 40 000) and HARPS 
data in standard HAM mode {R = 115 000). The spectra 
extraction was performed using the SOPHIE and HARPS 
pipelines. Following the techniques described by Baranne 
et al. (|T^M)) and Pepe et al. (I^UO^ . the radial veloci- 
ties were measured from a weighted cross-correlation of 
the spectra with a numerical mask. We used a standard 
G2 mask that includes more than 3500 lines. The result- 
ing cross-correlation functions were fitted by Gaussians to 
get the radial velocities and the associated photon-noise er- 
rors. The full width at half maximum of those Gaussians is 
12.6±0.2 kms"\ and its contrast is 27.7±0.8 % of the con- 
tinuum in the case of the HARPS data. The SOPHIE data 
gave similar parameters. We adjusted the number of spec- 
tral orders used in the cross-correlation in order to reduce 
the dispersion of the measurements. Indeed, some spectral 
domains are noisy, so using them degrades the accuracy of 
the radial velocity measurement. We finally used the orders 
10 to 38 for SOPHIE, and 5 to 71 for HARPS. 

Moonlight contamination was clearly visible in some 
spectra and in some cases at a radial velocity close to that 
of the target. Such contamination can affect the radial ve- 
locity measurements. Following the method described in 
Pollacco et al. (1^005]) and Hebrard et al. PM]) . we esti- 
mated and corrected the Moon contamination by using the 
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Fig. 4. Upper panel: Radial velocity measurements of 
CoRoT-18 with 1-a error bars as a function of time together 
with their Keplerian fit {top) and residuals of the fit. The 
fit is described below in Sect. |4.4| Lower panel: Same as 
above but as a function of the orbital phase. The data are 
from SOPHIE (red circles), HARPS (blue diamonds), and 
FIES (green squares). 



second optical-fiber aperture, which is targeted on the sky 
for both SOPHIE and HARPS, whereas the first aperture 
points toward the target. This induces radial velocity cor- 
rections up to 400 ms^^. 

The FIES observations were performed in January 2011 
under clear and stable weather conditions with seeing typi- 
cally in the range O'.'S— I'.'O. We used the 1'.'3 high-resolution 
fiber giving a resolving power R ~ 67000. Following the 
observing strategy described in Buchhave et al. (|2010p . 
three consecutive science exposures of 20 minutes were 
recorded for each observing night immediately preceded 
and followed by one long-exposed ThAr spectrum. Data re- 
duction and spectra extraction were performed using stan- 
dard IRAF routines. Finally, the FIES radial velocity mea- 
surements of CoRoT-18 were derived cross-correlating the 
science spectra with the spectrum of the radial velocity 
standard star HD 50692 (Udry et al. [TU^ observed with 
the same instrument set-up. 



The log of the observations and the radial velocity mea- 
surements are reported in Table [2] Radial velocity measure- 
ments have accuracies ranging between 23 and 98 ms~^ 
depending on the observing parameters. This table also 
shows the bisector spans that we measured on the cross- 
correlation functions in order to quantify the possible 
shape variations of the spectral lines. One SOPHIE spec- 
trum was too polluted by the Moon to allow any accurate 
bisector measurement. 

The radial velocity variations agree with Doppler shifts 
caused by a planetary companion, and the transit-signal 
detected from the CoRoT light curve could be interpreted 
as coming from a massive hot- Jupiter. We designate it as 
CoRoT-18b hereafter. 

The measurements are displayed in Fig. [4] together with 
their circular Keplerian fit, assuming the period and tran- 
sit epoch determined by the CoRoT light curve and refined 
with the photometric transit observed from the ground (see 
Sect. 3.4). The photometric and radial velocity data show 



good agreement. SOPHIE and HARPS radial velocities ob- 
tained with different stellar masks (FO or K5) produce vari- 
ations with the same amplitude as obtained with the G2 
mask, so there is no indication that their variations could 
be explained by blend scenarios implying stars of differ- 
ent spectral types. Similarly, the cross-correlation function 
bisector spans show neither variations nor trend as a func- 
tion of radial velocity (Fig. [s]). This reinforces the conclu- 
sion that the radial velocity variations are not caused by 
spectral line profile changes due to blends. 



• SOPHIE 

♦ HARPS 
■ FIES 





29.0 



29.5 
rodlol velocity (km/s) 



Fig. 5. Bisector span as a function of the radial velocities 
with 1-(T error bars. The ranges have the same extents in 
the X- and y-axes. The data are from SOPHIE (red circles), 
HARPS (blue diamonds), and FIES (green squares). 



3.3. Transit spectroscopy 

A transit of CoRoT-18b was observed in spectroscopy on 
2011 January 28. The goal was to detect the Rossiter- 
McLaughlin anomaly, which is an apparent distortion of 
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the stellar lines profile because of the transit of the planet 
in front of the rotating star. It allows the measurement 
of the sky-projected angle between the planetary orbital 
axis and the stellar rotation axis, usually denoted A (see, 
e.g., Bouchy et al. I2008p . The spectroscopic transit was ob- 
served with HARPS in the EGGS mode to improve the 
throughput. By co mparison to the HAM mode of HARPS 
used in Sect. 3.2 for the orbit determination (1" diame- 



ter fiber with a scrambler, allowing the resolution power 
R = 115 000 to be reached), the EGGS mode of HARPS 
uses a larger fiber (1'.'4) without scrambler. The spectral 
resolution is then reduced {R — 80 000) but the efficiency 
gain reaches a factor ^ 2 improvement. 

The target was continuously observed during a 5.5-hour 
sequence under good weather condition, with a seeing vary- 
ing between 0'.'7 and I'.'O. Twelve measurements were se- 
cured with exposure times ranging from 1200 to 1800 sec- 
onds, including five within the transit. The remaining ob- 
servations obtained before and after the transit are manda- 
tory for references. The radial velo cities were extracted as 
for HARPS/HAM data (Sect.[3]2|, but using fewer orders 
here for the cross-correlation (orders 18 to 69) to reduce 
the dispersion of the measurements. 

The HARPS/EGGS data are plotted in Fig. |6] The 
Rossiter-McLaughlin anomaly is detected, with an ampli- 
tude of ~ 100 ms~^, as expected according the rotation 
of the star and the depth of the photometric transit. The 
red shift during the first part of the transit and the blue 
shift during the second part indicate a prograde orbit. The 
symmetry of the feature agrees with an aligned system. 



3.4. Transit photometry 

A transit of CoRoT-18b was observed with the Euler 
1.2-m telescope on 2011 January 28 between 01:08 and 
05:36 UT, roughly at the same time and location as 
the Rossiter-McLaughlin ano maly observed in spectroscopy 
with HARPS (see Sect. 3.3 1. The goal was to refine the 
ephemeris. A total of 83 frames were recorded on the re- 
cently installed 4Kx4K E2V detector, with an exposure 
time of three minutes. As for the CoRoT light curve, we 
checked that this binning does not significantly affect the 
parameter retrieval by using the formalism presented by 
Kipping (j2010b|) . Standard calibration images were taken 
on the same night. The fluxes from the target and 20 refer- 
ence stars were extracted using standard aperture photom- 
etry with custom IDL routines. An average reference star 
was constructed by interactively selecting the stars that ex- 
hibited less real or instrumental variability. Nine stars were 
selected for this purpose. The final light curve was normal- 
ized to the median of the flux level after the egress of the 
transit. It is plotted in the lower panel of Fig.[3j 

We considered shot noise as a first estimation of the ac- 
curacy of the measurements, which was at the level of 10"'^ 
for this object. The accuracy level is later re-established 
by an evaluation of the reduced x^ of the fit process to 
take the corre lated noise present in the data into account 
(see Sect. 4.1 ). To better estimate the accuracy of the mea- 
surements, the dispersion of the data before and after the 
transit is 1.7 x 10"'^. The bump seen on the light curve 
near phase —0.01 could hint that the planet is transiting in 
front of a stellar spot; however, the amplitude of this event 
is within an the order of magnitude of the correlated noise 
so a spot detection cannot be claimed here. 




-0.02 0.02 

orbitol phase 



Fig. 6. Spectroscopic observation of the 2011 January 28 
transit of CoRoT-18b. Top: HARPS/EGGS radial veloc- 
ity measurements as a function of the orbital phase (filled 
diamonds), Keplerian fit ignoring the transit (dashed line), 
and final fit including the model of the Rossiter-McLaughlin 
anomaly (solid line). The vertical dotted lines show the 
times of mid-transit, first, second, third, and fourth con- 
tacts. Bottom: Residuals of the final fit. 



4. Analysis 

4.1. Light curves analysis 

4.1.1. Initial transit fit and light curve normalization 

To remove the modulations in the CoRoT light curve 
caused by rotating active regions on the stellar surface (see 
the top panel of Fig. [2]), we fit second-order polynomials 
to stretches of data on either side of each transit spanning 
approximately —2 to —1 and -1-1 to +2 times the transit du- 
ration. These sections were then normalized and kept for 
further light curve fitting, while flux measurements outside 
were discarded from the analysis. 

A detailed investigation of the immediate surroundings 
of the target (using Digital Sky Survey data) revealed that 
2.0 ± 0.1% of the flux in the photometric aperture was con- 
tributed by background stars. This was subtracted from the 
median out-of-transit flux before re-normalizing. The un- 
certainty in the contamination fraction translates into an 
effective uncertainty on the derived radius ratio due to con- 
tamination of ^ 0.0001. This is 15 times smaller than the 
flnal uncertainties derived below, so is negligible. 

The transits were then modeled using the formalism 
of Mandel & Agol p002p . with quadratic limb darkening 
coefficients Ui and U2 defined according to the standard 
law of the form /(/x)//(l) = 1 — mi(1 — ^u) — U2(l — /i)^, 
where /(I) and /(/i) are the specific intensities at the cen- 
ter of the stellar disk and at the angle 6 between the 
line of sight and the emergent intensity, respectively, and 
fjL = cos{9). We performed an initial least-squares fit (us- 



G. Hebrard et al.: CoRoT-18b, a massive hot jupiter on a prograde, nearly aligned orbit 



Table 2. Radial velocities of CoRoT-18. 



HJDuTC 


RV 


±lcr 


Bis. span 


exp. time 


S/N p. pix. 


Instrument/ 


-2 455 000 


(Jims"^) 


(kms^^) 


(kms'i) 


(sec) 


(at 550 nm) 


mode 


484.6704 


30.132 


0.041 


0.029 


3600 


16 


SOPHIE/HE 


485.6621 


29.030 


0.062 


0.184 


3600 


13 


SOPHIE/HE 


486.6469 


30.069 


0.051 


-0.196 


3600 


11 


SOPHIE/HE 


505.6367 


30.102 


0.032 


-0.074 


3600 


15 


SOPHIE/HE 


525.6027^ 


28.933 


0.086 


0.046 


2802 


8 


SOPHIE/HE 


527.5751''' 


29.006 


0.065 


-0.128 


3600 


11 


SOPHIE/HE 


529.61211^ 


29.178 


0.065 


0.117 


3600 


10 


SOPHIE/HE 


557.5227'' 


29.077 


0.096 


-0.209 


3600 


7 


SOPHIE/HE 


577.3831'' 


29.756 


0.041 


0.037 


3600 


16 


SOPHIE/HE 


578.3721'' 


29.241 


0.052 


0.165 


3600 


15 


SOPHIE/HE 


580.3673'' 


29.181 


0.094 


- 


3600 


13 


SOPHIE/HE 


538.7441 


28.979 


0.027 


-0.043 


3600 


11 


HARPS/HAM 


539.7253 


30.115 


0.038 


0.010 


3600 


9 


HARPS/HAM 


547.7838 


29.720 


0.023 


-0.011 


3600 


13 


HARPS/HAM 


548.6767'' 


29.324 


0.055 


-0.005 


3600 


7 


HARPS/HAM 


583.7515'' 


29.892 


0.051 


-0.072 


3200 


8 


HARPS/HAM 


590.7061 


29.732 


0.051 


-0.087 


3600 


7 


HARPS/HAM 


568.4512 


30.075 


0.066 


-0.049 


3600 


10 


FIES 


569.5086 


29.334 


0.070 


0.024 


3600 


8 


FIES 


570.5516 


29.881 


0.076 


0.031 


3600 


9 


FIES 
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f : measurements corrected from Moonlight pollution. 



ing the Levenberg-Marquardt algorithm) simultaneously to 
all 13 CoRoT transits, allowing the following parameters to 
vary: the period P, the epochs median of transit centers Tg, 
the planet-to-star radius ratio i?p/i?^, the impact param- 
eter b, the scaled semi-major axis a/ Ri,, and the combina- 
tions of the limb darkening coefficients Ui -I- U2 and wi — ^2- 
We considered a circular orbit. We discuss below the impact 
of possible low eccentricity on our results. The values for the 
limb darkening parameters obtained from this analysis were 
fully consistent with the values of ui = 0.47 and U2 = 0.21 
provided by Sing ([20T0|) for a star with Teg ^ 5500 K 
4.5, and [M/H 



logg 



-0.1 (see Sect. 



4.2 



remainder of our analysis, we fixed the limb dar 



). For the 
cening pa- 



rameters for the CoRoT light curve to these values. 



4.1.2. Detailed transit fit 

To estimate the photometric uncertainty of the CoRoT 
data, we calculated the standard deviation of the out-of- 
transit flux values, except those within 30 minutes of the 
ingress and egress of our initial fit. This provided an esti- 
mate of (7 = 1.7mmag, which was then used to perform 
a 100 000-step Metropolis-Hastings Markov Chain Monte 
Carlo (MCMC) analysis with the parameters P, Tq, &, 



a/ Ri,, and Rp/R^, allowed to vary, using our initial best-fit 
values from above as the starting points. We adjusted the 
jump scales for the free parameters until a step acceptance 
rate of 25-35% was achieved. The photometric uncertain- 
ties were then scaled up to cr = 1.9 mmag to give a reduced 
X^ of unity for the best-fit MCMC solution. This upwards 
scaling of the photometric uncertainties can be attributed 
to the presence of correlated noise in the CoRoT light curve 
(Pont et al. I2006p . often expressed as 



'N 



'W 



/N- 



2 



(1) 



where a\j is the variance between the mean flux values of 



'N 



bins containing N data points, and aw and cfr represent 
the "white" and correlated "red" components of the noise, 
respectively. The uncertainty on individual flux values is 
then given by cr = ajq for N = 1, so that a = \/(t| 



1-2 



The top panel of Fig. [7] shows calculated values for a as 
a function of TV for the CoRoT light curve. The variation 
between the mean flux values of bins containing only a sin- 
gle point {N = 1) is approximately 1.9 mmag, equal to the 
value that produces a reduced x^ of unity for the best-fit 
model. It also shows the expectation for a white-noise-only 
model that is forced to pass through the point at A^ = 10. 
For lower values of N , the white-noise-only model is not 
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CoRoT Inter-Bin RMS 




Euler Inter-Bin RMS 




Fig. 7. The standard deviation ctn of measured out-of- 
transit flux values as a function of the number of points 
A'' for the CoRoT light curve (top panel) and the Euler 
light curve (bottom panel). The dark curves indicate the 
best-fit models without red noise, i.e. only white noise. 



able to account for all of the variation in the data. This 
illustrates why it is necessary to scale the photometric un- 
certainties upwards from the standard deviation of the out- 
of-transit flux values. 

A similar procedure was then repeated for the Euler 
light curve. We first estimated the photometric uncer- 
tainty to be equal to the standard deviation of the out- 
of-transit measured fiux values, which was calculated as 
a = 1.9mmag. We then performed a 100 000-step MCMC 
analysis, allowing P, Tq, &, a/Ri,, and Rp/R-^, to vary. The 
limb darkening parameters ui and U2 were also allowed 
to vary, but this resulted in poor chain convergence. For 
this reason, we decided to set their values to Ui = 40 and 
1*2 = 30, as provided by Claret (2004) for a star with the 
same parameters as above. To investigate the effect of fixing 
the limb darkening parameters, we experimented with fix- 
ing them to other values provided by Claret ()2004p for stars 
with similar properties to CoRoT-18 and then performing 
a least squares fit to the light curve. In all cases, the fit- 
ted value for Rp/R^, changed, as expected, while the other 
parameters remained constant within the uncertainties. In 



particular, the fitted values for P and Tq were unaffected by 
varying the limb darkening parameters. This is important 
because the primary purpose of the Euler light curve is to 
refine the ephemeris of the orbit (see below) . We also ver- 
ified that the choice of limb darkening parameters for the 
Euler light curve did not affect the derived values for the 
other parameters in the simultaneous fitting to the CoRoT 
and Euler light curves, which is described below. This was 
done by ensuring that the fitted values for i?p/i?*, a/Ri, 
and b remained the same regardless of whether the Euler 
light curve was included in the fit. 

The bottom panel of Fig.[7]shows the variation between 
binned flux values as a function of bin number for the Euler 
light curve. It shows strongly correlated noise on time scales 
of ~ 10—20 minutes (between two to flve exposures ), which 
decreases significantly on longer timescales. This behavior 
is poorly modeled by Equation 1 above, with the solid line 
showing an illustrative fit for the white-noise-only case, i.e. 
(T/v = 0. Instead, we scaled the photometric uncertainties 
up to cr = 2.1mmag, which was the value required to give 
a reduced x^ of unity for the best-fit model. 

We next performed a joint MCMC analysis on the 
CoRoT and Euler light curves. This was done by initiating 
five independent chains at random locations in parameter 
space approximately ^ 5a away from the best-fit values 
determined from the initial fitting process. The period was 
held fixed to the value determined from the initial best- 
fit to the CoRoT light curve, but we allowed the transit 
midtimes to vary for each transit in order to inve stigat e 
the possibility of transit timing variations (see Sect. 4.1.3). 
Values for a/Ri,, Rp/Ri, and b were allowed to vary with 
the same values across all transits, and the limb darkening 
coefficients were held fixed to the values described above. 

The five chains were run in parallel for 200 000 steps. 
The initial 40 000 steps of each chain were then discarded 
to allow for the initial burn-in phase as the chains set- 
tled. We also checked when the x^ first dropped below the 
chain's median value, indicating the point at which the local 
neighborhood of the final solution had first been reached. 
Typically, this occurred within about 20 000 steps, imply- 
ing that the truncation we made at 40 000 steps is in fact 
conservative. The Gelman-Rubin statistic was then calcu- 
lated using the truncated chains for each of the free param- 
eters and, in all cases, found to be within approximately 1% 
of unity, indicating good mixing and convergence (Gelman 
& Rubin IT992P . The chains were then combined to obtain 
marginalized posterior distributions for each of the free pa- 
rameters, with medians and 1- a uncertainties reported in 
Table [3] and examples plotted in Fig. [8] The upper and 
lower Icr uncertainties respectively refer to the upper and 
lower bounds on the intervals containing 34.1% of the chain 
steps on either side of the median. 

A refined estimate for the orbital period was then ob- 
tained by fitting a model of the form 



Kin) ^Ta+: 



173 



(2) 



where TQ{n) denotes the transit midtime measured for the 
j^th Q].]-,iit^ Ji^ jg ii^Q epoch of the first CoRoT transit and 
Tb is the epoch of the Euler transit, which is 173 or- 
bital periods later. Then we ran two 200 000-step MCMC 
chains (one for Ta and one for Tb) and calculated the 
corresponding chain for the estimated period according to 
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Fig. 8. Examples of the results of the MCMC analysis with scatter plots showing the correlations between fitted param- 
eters. Values corresponding to the best-fit (lowest x^) solutions are indicated by red lines, while solid blue lines indicate 
the median values obtained for each parameter, and the dashed blue lines give the upper and lower 1-a uncertainties. 
The best and median solutions coincide for several parameters. 



F = (Tb - r4)/173. As expected from Pal (^003), the 
correlation between Ta and P chains was negative, while 
the correlation between the Tb and P chains was positive. 
We then calculated an optimal epoch using Equation 97 
of Pal (|2009p . which roughly corresponds to the median 
transit epoch, and find Topt = Ta + 32 x P. The new 
chain for this epoch led to the following median value 
for the optimal epoch, with associated 1-a uncertainties: 
To = Topt = 2 455 321.724f2±0.00018 HJD. By comparison 
to the classical computation using for example the epoch 
of the first transit, this method with the optimal epoch al- 
lows the uncertainty on Tq to be decreased by 2 sec, as well 
as the reduction of the correlation between Tq and P (see 
Fig. Is]) . The final ephemeris values are reported in Table Is] 
and plotted in Fig.[8J The best-fit solution for the combinea 
data set has a reduced x^ of 1.01. 

Finally, we estimated the effect of introducing a 
small eccentricity in the fit of the light curve. Indeed, 
Kipping (|2010ap has shown that assuming e = for an 
eccentric orbit could lead to underestimated uncertainties 
to a/R^, ratio and stellar density p*. We derive below the 
upper limit e < 0.08 at 95% confidence from the radial 



velocity measurements (see Sect. 4.4 1. We performed new 



tricity, with different possible values for the longitude uj of 
the periastron. Most extreme fits could provide a/R^ and 
/9* values at most to 2 a lower the values derived above. 
Thus we slightly decreased the final values on these pa- 
rameters and increased their uncertainties to account for a 
small possible eccentricity. 

4.1.3. Timing analysis 

The 13 successive transits observed by CoRoT and the 
observation of the additional transit with Euler offer 
an opportunity for transit timing variations (TTVs) re- 
search (Holman & Murray HUDSl). Figure [9] shows the TTVs 
measured by 



ATo{n) = T^{n)-To, 



(3) 



fits of the CoRoT light curve assuming this extreme eccen- 



where T'o(n) again denotes the transit midtime measured 
for the n"^ orbit and Tq denotes the refined transit midtime 
derived by fitting to the ephemeris equation, as described 
above. We obtain a reduced x^ of 1.3 for the hypothesis that 
there are no TTVs, i.e. Aro(n) = for n = 1, . . . , 14. When 
the ninth transit is removed (visible as the most discrepant 
point in Fig. ^ the reduce x^ improves to 0.98. Inspection 
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of this particular transit reveals that neither the ingress 
nor the egress was sampled, which could help explain why 
the measured TTV is somewhat discrepant. This is also 
the case for the fifth transit. Also, it is possible that the 
uncertainties in the measured transit midtimes are slightly 
underestimated because we did not account for unocculted 
spots that could result in different values for the effective 
radius ratio Rp/ R^, being measured from transit-to-transit 
as the spot coverage evolves; instead, we held Rp/Ri, fixed 
for all transits in the MCMC analysis. In either case, given 
the perfect agreement of the Euler ephemeris with all of 
the other CoRoT ephemerides, we conclude that there is 
no evidence of TTVs in the combined data set. Given that 
the characteristic uncertainty on our TTV measurements 
was about 60 seconds, we obtain a 3-cr TTV upper limit 
of 180 seconds. 



icant power at the first harmonic of the rotation period. 
We do not expect aliases here from the CoRoT satellite or- 
bital period, because it is well below the Nyquist sampling 
frequency of the data. We checked for that by binning the 
light curve to one point per orbital period (upper panel of 
Fig. [2]) to make sure that no systematics from the satellite 
rotation period can come through, and repeated the study. 
The results were identical. 

To refine the estimate of the rotation period, we fit a 
Gaussian to the periodogram around the main peak, shown 
in the lower panel of Fig. [2j and adopted the standard de- 
viation of the Gaussian as our period uncertainty, yielding 
^rot = 5.4 ± 0.4 days. This period is quite short for a late 
G-type star: even at the age of the Hyades, typical rotation 
periods for this spectral type are in the range of eight to 
nine days (Kawaler ,1989jl . 
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Fig. 9. Transit timing variation ATq for each transit of 
CoRoT-lSb. The 13 transits observed with CoRoT (on the 
left) were obtained about ten months before those observed 
with Euler (on the right). The 5"^ and 9*'' CoRoT transits 
do cover neither the ingress nor egress and therefore could 
conduct to erroneous timing measurements. 



4.1.4. Rotation period 

The CoRoT white-light curve of CoRoT-18 shown in the 
top panel of Fig.[2]displays obvious signs of rotational mod- 
ulation of star spots, with at least two large active regions 
present on the stellar surface at any one time. To estimate 
the rotation period we first cut out a section lasting twice 
the transit duration around each transit, then fit a straight 
line to the out-of-transit light curve to remove any long- 
term trend. We also corrected for a discontinuity near the 
end of the light curve by subtracting ISmmag from all data 
taken after the discontinuity. This detrending procedure is 



slightly different from the one presented above in Sect. 4.1 



We estimated the rotation period by fitting sinusoids 
at 1000 trial periods ranging from 1 to 15 days to the cor- 
rected out-of-transit data. The resulting best-fit amplitudes 
are shown as a function of trial period in the bottom panel 
of Fig. [2j The dominant peak clearly corresponds to the 
interval between repeat appearances of individual active 
regions, hence to the rotation period. There is also signif- 



4.1.5. Planetary occultation upper limit 

We derived the depth upper limit of a possible occultation 
of the planet by the star by fitting a Mandel & Agol (|2002|) 
transit model at phase 0.5, with no limb darkening and a 
transit depth reduced by a factor Fp/F^,, where Fp is the 
planetary flux and F^ the stellar flux. All parameters were 
kept fixed to the best-fit values derived from the transit, 
except for the eclipse depth, for which we obtain a best-fit 
value of 0.05 ± 0.19 mmag, i.e. consistent with zero. The 
resulting 3-cr upper limit is 0.61 mmag. This is not partic- 
ularly stringent, as one can expect a planetary to star flux 
ratio of 0.40 mmag in the optimistic case of a geometric 
albedo equal to one. 



4.2. Stellar analysis and classification 

To determine the planetary parameters with an as high as 
possible precision we need to know the physical conditions 
of its host star. Seager & Mallen-Ornelas (|2003p . among 
others, have shown that ideally there is one stellar parame- 
ter, the stellar density, which can be obtained from a tran- 
sit light curve of sufficient photometric precision. From this 
parameter, it is possible (with a number of assumptions) 
to derive, through modeling, other physical parameters of 
the system. Nevertheless, as pointed out in this context by 
e.g. Fridlund et al. (2010), high-precision photometric and 
spectroscopic measurements that has been carried out on 
other exoplanet host stars, do suggest that this rarely infers 
reliably to the other properties of the star - mainly because 
of flaws in stellar theory (Winn et al. I2008P . 

We used the two sets of HARPS observations to perform 
this analysis: the HARPS/HAM data (co-additi on o f six 
spectra totaling 5.8 hours of integration, see Sect 3.2 ) and 
the HARPS/EGGS data (co-addition of 12 spectra total- 
ing 5.4 hours of integration, see Sect 3.2). Due to its lower 
resolution power (80 000 vs. 115 000)~the HARPS/EGGS 
data presents the higher signal-to-noise ratio (~ 35 in the 
continuum at Ha) . We analyzed both sets of observations 
and find no significant differences in the stellar parameters 
beyond the internal 1-ct error. Three observations are im- 
mediately made while inspecting the co-added spectra. The 
appearance is that of a cool star, the line profiles are rela- 
tively broad (i;sini*= 8.0±1.0 kms""'^), and there is a faint 
absorption (equivalent width of ~ 40 mA) at the location 
of the Li I (6707.8 A) line. There is no obvious detection of 
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any Can chromospheric emission. This wsinz,t direct mea- 
surement agrees with this derived from the cross-correlation 
function following the Santos et al. (|2002p methodology. 

To determine the spectroscopic parameters, we used 
the Spectroscopy Made Easy (SME, version 162, February 
2011) software package (Valenti & Piskunov 119961 Valenti 
& Fischer |2005() . SME calculates synthetic spectra and fits 
the observations to it. All the normal stellar parameters 
(Teff, \ogg, wsini*, [Fe/H], abundances, etc.) can be used 
either as input or as free parameters to solve for. A grid 
of stellar models (Kurucz models) was utilized in order to 
determine the fundamental stellar parameters iteratively. 
This was achieved by fitting the observed spectrum directly 
to the synthesized spectrum and minimizing the discrep- 
ancies using a nonlinear least-squares algorithm. SME re- 
quires atomic line data in order to synthesize a spectrum. 
We utilized input from the Vienna Atomic Line Database 
(Kupka et ah . 1999. Piskunov et al.\TMB- 

Using SME and a sample of more than 1000 stars, 
Valenti & Fischer (|2005p found uncertainties of 44K in Toff, 
0.06 dex in \ogg, and 0.03 dex in [M/H] per measurement. 
Based on the CoRoT material (stars hosting CoRoT plan- 
ets or candidates that have not been positively shown to 
host planets), we find slightly larger errors than Valenti & 
Fischer (|2005p : 70 - 100 K in Toff, depending on the signal- 
to-noise in the continuum of the spectrum at the location 
of the Balmer lines, 0.05 — 0.1 dex in \ogg, depending also 
on the spectral type and on which ion we used, and fi- 
nally 0.05 — 0.1 dex in [M/H]. However, by comparing the 
measurements with model isochrones they found a larger, 
systematic offset of ~ 0.1 dex and a scatter that can oc- 
casionally reach 0.3 dex in logs'. In. CoRoT-18, we find an 
internal discrepancy using SME of 0.1 dex. We therefore 
assign 0.1 dex as our l-cr precision. 

We found Teff = 5U3K ± 100 K from the profile of the 
Balmer lines. We determined the metallicity and found con- 
sistent results from different ions indicating a star of slightly 
lower than Solar metalhcity: [M/H] ~ -0.1 ±0.1. The log^ 
was determined utilizing the lines of Cai, Mgi and Nai, 
finding a consistent result of 4.4 ±0.1. The spectroscopic 
parameters of CoRoT-18 are summarized in Table [3] 



evolution (see Figs. 2-4 from Guillot & Havel l20lT|l . We 
used the CESAM evolution code (Morel & Lebreton I^DU^ 
to calculate these evolution tracks. In Fig. 10 we plot the 
solutions in the stellar-mass, age parameters, with col- 
ors that depend on their quadratic distance to the effec- 
tive temperature-stellar density constraints. These point 
towards either a young age, less than 30 Ma (at 3a), or 
an old one, more than 8 Ga at la and more than 4 Ga at 
3(7. The situation is thus different than for CoRoT-2, for 
which a continuum of young and late-type solutions was 
found (see Guillot & Havel [MDl- For CoRoT-18, these so- 
lutions contradict with the other age indicators. 
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Fig. 10. Constraints obtained from stellar evolution models 
for the age and mass of CoRoT-18. The colored circles cor- 
respond to constraints derived from stellar evolution mod- 
els matching the stellar density and effective temperature 
within a certain number of standard deviations: less than 
1(7 (green), 2a (blue), or 3ct (yellow). 



4.3. Stellar evolution tracks and the age problem 

Altogether, CoRoT-18 seems remarkably similar to CoRoT- 
2 (Alonso et al. [MTH]) : the two stars (CoRoT-18 vs. 
CoRoT-2) have comparable effective temperatures (5440 
vs. 5450 K), metallicities (—0.1 vs. 0.0), spin periods (5.4 
vs. 4.5 days), and vsini^, (8.0 vs. 12 kms~^), and they are 
both active, with peak-to-peak photometric variabilities of 
~ 2% and ^ 4%, respectively. In addition CoRoT-18b and 
CoRoT-2b are the only known planets (transiting or not) 
orbiting a star colder than 6000 K and with a large wsini*. 
The other planet-host stars in this temperature range all 
have iisini* values in the range [0 — 5] kms^^. All the 
planet-host stars having wsini* ^ 10 kms~^ or larger are 
F-stars - except CoRoT-18 and CoRoT-2. 

However, the inferred stellar densities for CoRoT-18 and 
Co RoT-2 (1.35±0.25 vs. 1.814+°;[]|°gcm-3 - see Gillon et 
al. I2010p slightly differ. As a result, the effective tempera- 
ture, metallicity, and density constraint for CoRoT-18 are 
consistent with evolution tracks for solar-mass stars that 
are either particularly young and still on the pre-main se- 
quence or old and towards the end of the main-sequence 



First, as for CoRoT-2, the rapid rotation of the star 
points towards a young age. Bouvier (2007) has compiled 
hundreds of rotational period measurements from photo- 
metric surveys of young open clusters with ages up to 
625 Ma (Hyades), and used them to model the rotational 
evolution of stars in several mass bins. A range of models 
is needed to reproduce the data at any given age, and con- 
straints are scarce beyond 200 Ma for stars with masses be- 
tween 0.8 and IMq. However, extrapolation of the models 
calibrated at earlier ages suggests that such stars are not ex- 
pected to retain rates ~ 5 times faster than the Sun (which 
is the case for CoRoT-18) beyond ages of 500 -600 Ma. 

Second, the lithium equivalent width favors the latter 
end of this age range. According to a recent compilation of 
Li depletion measurements, also from stars in young open 
clusters (Hillenbrand et al. I2009p . the equivalent width of 
40 mA measured for CoRoT-18 is typical of stars of this 
spectral type at the age of Ursa Majoris (500 Ma) or the 
Hyades (625 Ma), while typical equivalent widths for similar 
stars in M34 and M7 (- 200 Ma) are about 100 mA. This 
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would thus instead suggest an age of several hundred Ma 
for CoRoT-18. 

Thus, for now we are unable to estimate the age of 
CoRoT-18, even if it seems to be young. This illustrates 
the difficulty in determining the age of stars. We adopt 
the stellar mass M^, = 0.95 ± 0.15 M0. The conservative 
error bar is large enough to agree with both pre-main- 
sequence case and old star main-sequence star. This im- 
plies a semi-major axis a = 0.0295 ±0.0016, a stellar radius 
R^, = 1.00 ± 0.13 R0, and thus a planetary radius Rp — 
1.31±0.18Rjup. The equilibrium temperature of the planet 
assuming an isotropic zero-albedo i s Teg = 1550±90K. The 
rotation period 5.4±0.4 days (Sect. 4.1.4) is consistent with 
the high value vsini^, = 8.0 ± 1.0 kms~^ and is fast for a 
late G-type star. Using i?^, = 1.00 ± 0.13 R©, this trans- 
lates into an inclination of the stellar rotation axis that is 
i* = 70° ± 20°, so the star is seen nearly edge-on. 

4.4. Radial velocities fit 

We fitted the radial velocities with a Keplerian model. The 
period and the epoch of the transit were fixed to the val- 
ues obtained from the light curves analysis (Sect. 4.1). If 
the relative accuracy of the radial velocity measurements is 
high (a few tens of ms~^ here in the case of CoRoT-18), 
their absolute accuracy in heliocentric or barycentric frames 
could be ten times less good, so a radial velocity shift was 
free to vary in the fit between the three datasets used for 
the orbit (SOPHIE, HARPS/HAM, and FIES), and we fi- 
nally obtained one systemic radial velocity for each of the 
three instruments. The HARPS/EGGS data secured during 
a transit do not significantly constrain the orbit, and they 
are used below for the analysis of the Rossiter-McLaughlin 
anomaly. We found the upper limit e < 0.08 at 95% confi- 
dence for the eccentricity of the orbit and thus assumed a 
circular orbit, as usually is the case for hot jupiters. In the 
case of a slightly eccentric orbit, its orientation is not well 
constrained, with the longitude of the periastron included 
in the range -60° <uj < 100°. 

The final fit of the radial velocities is plotted in Fig. 4 
The derived orbital parameters are reported in Table 3 
together with error bars that were computed from x^ vari- 
ations and Monte Carlo experiments. The radial velocity 
variations present a semi-amplitude K = 590 ± 14 ms~^, 
corresponding to a planet with a mass Mp — 3.47 ± 
0.38 Mjup. This assumes M^ = 0.95 ± 0.15 Mq for the host 
star, which here is the main source of uncertainty on Mp. 

The standard deviation of the residuals to the fit is 
ao-c = 41.0 ms~^ for the whole dataset (35.9, 44.4, and 
43.0 ms-i for SOPHIE, HARPS/HAM, and FIES, respec- 
tively). The reduced x^ is 1-02 for the 22 radial veloci- 
ties used in the fit. We do not detect any drift over the 
106-day span of the radial velocity, with an upper limit of 
200 ms~^ a~^ at 95% confidence. We can thus exclude any 
extra planet in the system with a mass higher than 3 Mj^p 
and a period shorter than 200 days. 



4.5. Planetary evolution 

We have seen that CoRoT-18 and CoRoT-2 arc similar for 
what concerns their stars. Their planets (CoRoT-18b vs. 
CoRoT-2b) are also strikingly similar, in terms of orbital 
periods (1.90 vs. 1.74 days), masses (3.4 vs. 3.7Mjup), and 



equilibrium temperatures (1550 vs. 1539 K). The inferred 
planetary radii are 1.31 ± 0.18 Rjup for CoRoT-18b and 
1.466tg;o44Rjup for CoRoT-2b (Gillon et al. lMTOl) . 
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Fig. 11. Constraints obtained on the age and radius of 
CoRoT-18b. The colored circles correspond to 1 a (red), 2 a 
(yellow), or 3(7 (blue) solutions. The figure also plots the 
same constraints for CoRoT-2b (grayscale), when not in- 
cluding the effects of spots (see Czesla et al. 2009j Guillot & 
Havel ECTTTj) . The evolution tracks show the progressive con- 
traction of a 3.5Mjup planet with Tcq ~ 1600 K, in the so- 
called "standard approach" , when increasing atmospheric 
opacities by a factor 3.5, and when dissipating 10^^ ergs~^ 
at the center (see Guillot & Havel l20Tl1 for a description of 
the models). 



Figure [TT] compares the planetary radii obtained as a 
function of age for the two planets, using the approach de- 
scribed in Guillot & Havel pOlip . To use the same ap- 
proach for the two planets, when calculating the size of 
CoRoT-2b, we did not account for the effect of spots (see 
Czesla et al. l2009l Guillot & Havel [20TT]) - in the case 
of CoRoT-2b, these yield a ~ 3% increase in the inferred 
radius. For CoRoT-2b, the age constraints from evolution 
models matching (Tcff, /?*, \ogg, [Fe/H]) are weak and 3tT 
solutions are found anywhere between 20 Ma to more than 
15 Ga). The inferred radius is extremely large for an object 
of this mass and planetary evolution models require a re- 
cent (~ 20 Ma) dramatic event (birth, giant impact, close 
encounter, and circularization of the orbit) to explain it. 
Because of its lower stellar density, solutions for CoRoT-18 
are found at both extremes of the age range (10 to 25 Ma 
and more than 4Ga, see Sect. 4.3). The constraints on the 
size of CoRoT-18 are weaker so that they cannot be used 
to confirm or refute possible inflation mechanisms or pos- 
sible compositions, at least as long as the cause of the age 
mismatch betwe en st ellar evolution models and youth indi- 
cators (see Sect. 4.3 ) is not found. 



If the mismatch found for CoRoT-18 between the var- 
ious age indicators can be linked to the poorly known 
physics of young stars, it may lead to completely reconsider 
the problem of the inflated size of CoRoT-2b and other ex- 
oplanets orbiting young stars. 
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4.6. Rossiter-McLaughlln anomaly analysis 

The radial velocity measured during the 2011 January 28 
transit were fitted in order to derive the sky-projected angle 
A between the planetary orbital axis and the stellar rota- 
tion axis. We applied the analytical approach developed 
by Ohta et al. (12005]) to model the Rossiter-McLaughlin 
anomaly shape, which use here ten parameters: the stel- 
lar limb-darkening linear coefficient e, the transit param- 
eters Rp/Ri,, a/R^, and i, the parameters of the circular 
orbit (P, To, and K), the HARPS/EGGS systemic radial 
velocity, and finally vsini^ and A. We adopted e — 0.722 
computed by Claret ()2004p in the g' filter corresponding 
to the HARPS wavelength range. The transit and orbital 
parameters were determined above from the light curves 
and radial velocity fits, and their uncertainties are negligi- 
ble for the fit of the Rossiter-McLaughlin anomaly shape, 
according the uncertainties of the HARPS/EGGS radial ve- 
locities. The main parameters that play a role in this fit are 
the systemic velocity. A, and wsini*. As these parameters 
are correlated in the Rossiter-McLaughlin fit, we computed 
the x^ of the fit on a three-dimensional grid scanning their 
possible values. 

The systemic velocity we obtained for the transit ob- 
served with HARPS/EGGS, 29.550 ± 0.016, is similar to 
this obtained in Sect. 14.41 for the orbit observed with 
HARPS/HAM. It is constrained thanks to the observations 
secured immediately before and after the transit. The con- 
fidence interval contours estimated from x^ variations for 
the A and wsinz* are plotted in Fig. [12] We thus obtained 
A= -10° ±20° and w sin i* = 8. 5 ±2. 5 km s~^ The best fit 
is plotted in Fig. [6] The dispersion of the residuals to the 
fit is 46.5 ms~^. This agrees with the expected error bars 
on the radial velocities, and this is similar to the dispersion 
of the residuals to the fit of the orbit presented in Sect. 4.4 




-60 -40 -20 20 40 

A (degrees) 

Fig. 12. x^ isocontours for our modeling of the Rossiter- 
McLaughlin effect as a function of A and wsini*. The dia- 
mond shows the lowest-^^ value. 



The Rossiter-McLaughlin anomaly detection only 
stands on a few points; however, the shape of the ra- 
dial velocity variations during the transit agree with a 
Rossiter-McLaughlin feature, with A ~ 0° and the ex- 



pected amplitude. As a statistical test for the Rossiter- 
McLaughlin anomaly detection, we computed the x^ over 
the 12 measurements secured during the transit night, and 
we get 12.5 and 24.7 for the fits including or not the 
Rossiter-McLaughlin anomaly, respectively. Including the 
Rossiter-McLaughlin model in the fit thus implies a fac- 
tor two improvement in the x^, for basically two extra 
free parameters (A and wsini*, which mainly constraints 
the Rossiter-McLaughlin shape). We performed an F-test 
which indicates there is a probability > 70% that the x^ 
improvement actually is due to the Rossiter-McLaughlin 
anomaly detection. 

Usually only the sky-projected value A of the obliq- 
uity ip could be measured because the inclination i* of 
the stellar rotation axis rem ains unknown. Here we esti- 
mated i* = 70° ± 20° (Sect. [42]), so the sky-projected A 
could be translated into the actual obliquity. We obtain 
ip = 20° ± 20°. This value remains inaccurate, due to the 
significant uncertainties on A and z* . It allows, however, the 
conclusion that the orbit of CoRoT-18b is prograde and 
nearly aligned. This additionnaly reinforces the similarity 
between the CoRoT-18 and GoRoT-2 planetary systems, 
since CoRoT-2b is also prograde and aligned (Bouchy et 
al.EPOS). 

The V sini* value obtained from this fit agrees with this 
obtained in Sect. |4.2| from spectral analysis, usini* = 8.0 ± 
1.0 kms~^. As discussed by, e.g., Hirano et al. POIOP and 
Simpson et al. (|201ip . modeling the Rossiter-McLaughlin 
anomaly could produce biased v sin i* measurements when 
rotational broadening of the stellar lines is significantly 
larger than the instrument resolution. We do not see that 
effect here, possibly because of the long exposure times of 
each exposure. 

5. Conclusion 

We reported the detection of the 18"' transiting exoplanet 
detected by the CoRoT project. This giant planet was 
discovered thanks to the high-accuracy, continuous pho- 
tometry obtained by the CoRoT satellite and the photo- 
metric and spectroscopic follow-up performed on ground- 
based telescopes. CoRoT-18b is a massive hot jupiter or- 
biting a faint G9V star. Its mass is Mp — 3.47 ± 0.38 Mjup, 
and its radius Rp — 1.31 ± 0.18 Rjup, implying a density 
Pp = 2.2 ± 0.8g/cm'^. The period of the circular orbit is 
1.9000693 ± 0.0000028 days. It is known with an accuracy 
better than 0.25 seconds and the mean mid-transit epoch 
with an accuracy of 20 seconds. The mass of the host star 
is A'U = 0.95±0.15Mq and its radius R^ = l.OO±O.13R0. 
The parameters of this system are summarized in Table [3] 
The parameters of CoRoT-18b are similar to those of 
CoRoT-2b (AlonsoEOOHl Gillon et al.[2 010]), a nd to a lesser 
extent to those of Co RoT-l lb (Gandolfi [50101) and CoRoT- 
17b (Csizmadia et al. I201ip . the other massive hot jupiters 
found with CoRoT. Interestingly CoRoT-18b is found to be 
either particularly young (a few tens of Ma) or old (> 4 Ga) 
from stellar evolution models matching the star's effective 
temperature and inferred density, but according both to 
its lithium abundance and to its relatively fast rotation, 
it would be expected to be modestly young (several hun- 
dred Ma). This mismatch potentially points to a problem 
in our understanding of the evolution of young stars, with 
possibly significant implications for stellar physics and the 
interpretation of inferred sizes of exoplanets around young 
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Table 3. Planet and star parameters. 
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Fitted transit parameters 


Planet orbital period P [days] 


1.9000693 ± 0.0000028 


Transit center epoch To [HJD] 


2 455 321.72412 ±0.00018 


Scaled semi-major axis a/Rt 


6.35 ± 0.40 


Radius ratio k = Rp/R, 


0.1341 ±0.0019 


Impact parameter b = a cos i/Rt 


oAotrA 


Deduced transit parameters 




Orbit inclination i [°] 


86.5li:^ 


Transit duration Ti_4 [h] 


2.387 ±0.037 


Egress/ingress duration T1-2 ~ 73-4 [h] 


0.331 ±0.038 


mI^^ /Rt [Solar units] 


0.99 ±0.06 


Stellar density p» [g cm~ ] 


1.35 ±0.25 


Results from radial velocity observations 




Radial velocity semi- amplitude K [ms~^] 


590 ± 14 


Orbital eccentricity e 


< 0.08 


SOPHIE systemic velocity [kms~^] 


29.533 ±0.016 


HARPS/HAM systemic velocity [kms^^] 


29.572 ±0.015 


FIES systemic velocity [kms^^] 


29.657 ±0.034 


O-C residuals [ms~^] 


41.0 


Spectroscopic stellar parameters 




Effective temperature Teff [K] 


5440 ± 100 


Surface gravity logg [dex] 


4.4 ±0.1 


Metallicity [Fe/H] [dex] 


-0.1 ±0.1 


Rotational velocity w sin u [kms~^] 


8.0 ±1.0 


Spectral type 


G9V 


Star mass [Mq] 


0.95 ±0.15 


Stellar and planetary parameters 




Star radius [Rq] 


1.00 ±0.13 


Distance of the system [pc] 


870 ± 90 


Stellar rotation period Prot [days] 


5.4 ±0.4 


Stellar inclination i* [°] 


70 ±20 


Orbital semi-major axis a [AU] 


0.0295 ±0.0016 


Planet mass Mp [Mj ] 


3.47 ±0.38 


Planet radius -Rp[Rj] 


1.31 ±0.18 


Planet density pp [g cm~'^] 


2.2 ±0.8 


Equilibrium temperature Toq [K] 


1550 ± 90 


Rossiter-McLaughlin parameters 




Sky-projected obliquity A ["] 


10 ±20 


HARPS/EGGS systemic velocity [kms"^] 


29.550 ±0.016 


O-C residuals (HARPS/EGGS) [ms^i] 


46.5 


Obliquity -ip [°] 


20 ±20 



stars. In addition, CoRoT-18b and CoRoT-2b are the only 
known planets (transiting or not) orbiting a fast-rotating, 
G-type or cooler star. 

The orbit of CoRoT-18b is prograde, with a spin-orbit 
angle V' = 20° ± 20° (sky-projected value A = -10° ±20°) 
hence an obliquity in agreement with 0. Schlaufman (|2010p 
and Winn et al. (j2010b|) have shown that misaligned plan- 
ets tend to orbit hot stars. With an effective temperature of 
5440 ± 100 K for the host star, the CoRoT-18 system sup- 
ports this trend. Exceptions to this trend are increasing, 
however (Moutou et al. lMTTl Brown et al. [20TT|) . 

Hebrard et al. (2010j 2011) have hypothesized that most 
of the massive planets are prograde and moderately but sig- 
nificantly misaligned, whereas the less massive planets are 



distributed in two thirds of the prograde, aligned systems 
and one third of the strongly misaligned systems. Being pro- 
grade and nearly aligned, CoRoT-lSb is at the upper limit 
of the low-mass range. Interestingly, a limit is also appar- 
ent in the mass distribution of known transiting planets. 
This is shown in Fig. |13[ where one can see a decreasing 
abundance of planets with increasing planetary mass up 
to 4.5Mjup. No transiting planets are known in the range 
Mp= [4.5 — 7] Mjup, and a few are known with Mp> 7 Mjup. 
It is difficult to imagine that a bias is provoking the lack 
of massive planets, as they are easier to detect. This differ- 
ent distribution suggests that planets below 4.5 Mjup could 
have a different nature or history than those above 7Mjiip. 
This is reinforced by the different obliquity distribution. 
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Fig. 13. Histogram of the number of known transit- 
ing planets as a function of their mass (data from 
littp://exoplanet.eu ). 
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